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ABSTRACT 


The Navy’s primary analysis of damage control and stability to date has been 
under static conditions. Dynamic effects, such as progressive flooding, and the dynamic 
damage control procedures, such as hole patching and dewatering, have not been 
included in present design requirements. The goal of this thesis is to develop and test a 
stand-alone progressive flooding model. This model can be used to evaluate the transient 
and steady state characteristics of shipboard progressive flooding. Several improvements 
Over previous studies are introduced and their effects are assessed. A sensitivity analysis 
study is performed through a systematic series of runs for a variety of hull forms. These 
results can be used to aid engineers of future ship designs in the use of damage control 


techniques and parameters. 
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I. INTRODUCTION 


an BACKGROUND 


Damage control practices are one of the most effective ways to increase ship 
survivability. Damage control practices can be categorized in two ways: static, such as 
the use of watertight bulkheads and other watertight fixtures and dynamic, such as the use 
of repair and dewatering techniques. Static damage control practices have been used for 
centuries as early shipbuilders designed watertight bulkheads into their design. Although 
it can not be documented it is suspected that the use of dynamic damage control 
procedures dates back even further since common sense dictates the use of dewatering 
and repair techniques on a damaged vessel. In the 1930’s the U.S. Navy began, for the 
first time, implementing damage control and damage stability criterion into future ship 
design. After World War I the Navy began conducting damaged stability studies on new 
combatants. These studies led to some of the damage control practices implemented 
during World War II [Ref. 1]. Post World War II studies conducted by the Bureau of 
Ships, the current Naval Sea Systems Command (NAVSEA), resulted in the requirement 
that Naval vessels be able to withstand a hole that was equal to 15% of there length for 


combatants and 12.5% of there length for auxiliaries [Ref. 2]. 


B. CURRENT DESIGN REQUIREMENTS AND PROCEDURES 


The current standards are delineated in NAVSEA Design Data Sheet 097-1 (DDS 
097-1). These standards are based upon a study by T. H. Sarchin and L. L. Goldberg and 
the BuShip study of 1947. The DDS 097-1 criterion for compartmentation of category I 
ships (which includes combatants) is that the ship withstand a rapid influx of water from 
a shell opening equal to 15% of the ship’s length at any point along the length of the ship. 
Reserve buoyancy requirements are that the equilibrium line not be above the margin 
line, which lies three inches below the main deck [Ref. 1]. The Navy currently uses a 
naval architecture program called the Ship Hull Characteristics Program (SHCP) as its 
primary tool to implement these requirements [Ref. 3]. This is a Fortran based program 
that uses a geometry interpreter and several naval architecture functions to evaluate the 
hull-form. This procedure is limited because it does not implement any dynamic damage 
control techniques such as dewatering or repair. This evaluation also does not take into 


account the effects of progressive flooding into adjacent compartments. 


GC: THESIS OUTLINE AND CONTRIBUTIONS 


By the year 2010, the Chief of Naval Operations (CNO) has endorsed a series of 
operational characteristics that must be incorporated into surface combatants. One of the 
primary characteristics is that the ship must retain the ability to conduct combat 
operations, even though it has sustained hull damage and flooding [Ref. 1]. This 


requirement has further outdated the static damage control and stability tests that are 


conducted with the SHCP. David Taylor Research Center has conducted dynamic 
damage control tests on certain surface combatants to evaluate their stability. Further 
dynamic testing of surface combatants 1s needed as the Navy shifts to more performance 
based requirements. This thesis designs a progressive flooding tool that will aid to 
further implement some of the operational characteristics of surface combatants. Some 
of the advantages and limitations of this progressive flooding simulation program are: 

1. Visual Basic allows the development of a portable analysis tool that can be taken 
to any location that has a portable computer with a Windows 95/98/NT operating 
system. The simulation model can easily be installed on a portable computer in a 
matter of minutes and no supporting software 1s required. 

2. Visual Basic allows the programmer of the analysis tool to include many of the 
graphical and visual features that further analyze transient and steady state 
characteristics of the progressive flooding process. Furthermore, these features 
are implemented in a familiar and user-friendly interface, which will keep the user 
learning curve at a minimum. 

3. The user can conduct real-time damage control techniques and include them into 
the simulation, such as hole plugging and dewatering with installed pumps. 

4. The user can decide what will be the limiting criterion for the progressive 
flooding simulation such as available freeboard and minimum transverse 
metacentric height. 

5. The Visual Basic language allows the user to develop an accurate model, which 
includes all necessary naval architecture parameters. Extensions to include 


dynamics effects such as water sloshing are possible and can be easily 


incorporated into the program. On the other hand, it will be rather cumbersome to 

include an accurate simulation of pump characteristics and the ship’s fire main 

system. If such studies are needed, we recommend the use of a SIMSMART 
based tool [Ref. 1] that includes all necessary algorithms for solving a network of 

Bernoulli and continuity equations. 

Chapter II of this Thesis presents the mathematical model that was utilized in the 
program. The basic flow-rate equations along with the fundamental naval 
architectural equations are presented for a rectangular barge, a Wigley hull, which is 
an analytical hull form explained in Chapter II, Section B, and an arbitrary hull form. 
The development and operation of the main simulation program are also described in 
this chapter. Chapter III presents results on program validation. This is accomplished 
by comparing the results for the Wigley hull with results obtained by treating the 
same hull as if it were an arbitrary hull with user-provided hydrostatic data. Based on 
the results of various simulation scenarios, several trends are identified and discussed. 
Finally, in Chapters IV and V, we summarize the main conclusions from this study 


and offer recommendations for further research. 


Il. DEVELOPMENT OF THE SIMULATION MODEL 


A. APPROACH OF MODEL DEVELOPEMNT 


The primary reason for development of this model is to analyze the transient and 

steady state characteristics of surface ship progressive flooding. When a hull form 
composed of various compartments suffers damage and becomes open to the sea the 
compartment that is open to the sea floods. If the damage was not severe and the 
watertight integrity remains intact then the flooding is isolated to only that compartment. 
However, it is more realistic to consider the damage severe enough to disrupt the 
watertight integrity of that compartment and introduce progressive flooding into the 
surrounding compartments through small holes or fragmentation. Progressive flooding 
may also occur due to improper maintenance or normal wear and tear of watertight 
fixtures. It can eventually become serious enough to cause the ship to founder (sinking 
caused when the remaining buoyancy is less than the ship’s weight) or the ship can 
capsize due to loss of stability. 

The main focus of this model will be on three compartments. The user will define 
the bulkhead locations at the beginning of the simulation. The primary compartment will 
contain the hole(s), which open it to the sea and allow it to flood. This compartment can 
have as many as two holes open to the sea. The user will define the diameter, height and 
flow coefficient of al! holes at the beginning of the simulation. Default values for the 
flow coefficients (equal to 0.7) are provided by the program. The two additional 


compartments will be called secondary compartments. They are located immediately 


forward and aft of the primary compartment. The secondary compartments can contain 
as many as three secondary fragmentation holes, which allow the primary compartment 
to progressively flood into the secondary compartments. 

The model for this analysis was developed using the Microsoft Visual Basic 
programming language. One of the main goals for this thesis was to develop an analysis 
tool that was independent of any other programs. It was also desirable for the analysis 
program to be easily exportable so that it may run on other platforms. Visual Basic is an 
excellent tool for developing stand-alone, exportable programs. Once a program is 
developed it can be compiled into an exportable package that can be easily installed on 
any desktop computer. Visual Basic also offers the programmer access to some of the 
mouse driven visual features common in most modern programs, allowing the end user of 
the analysis program to demonstrate a real-time, interactive scenario. Some of the real- 
time features, which are included in this simulation tool, are hole clogging (or repair) and 
additional pumping to aid in the dewatering process. 


The user can use this program to monitor various aspects of shipboard progressive 
flooding. Some of the parameters of interest are metacentric height (GM ), available 


freeboard, and trim forward and aft. The GM is a measure of the ship’s transverse 


stability and can be continuously observed while the simulation is running. A significant 


decrease in GM may allow the ship to capsize when a list is encountered. The available 
freeboard can also be continuously monitored while the simulation is running. The 
available freeboard is a measure of the reserve buoyancy that the ship has at an instant in 
the simulation. When the freeboard is depleted the ship will founder. The trim that the 


ship experiences during the simulation may be of interest to the user if the hull has certain 


trim restrictions. The simulation will automatically terminate when the simulation 


reaches a user defined minimum GM orall available freeboard has been depleted. All 
simulations in this version are performed using an explicit Euler’s formula, which is 
adequate for the time scales considered. In case more elaborate effects are required, such 
as compartment water sloshing, the method can be easily changed into a Runge-Kutta 


based formula. 


B. HULL FORMS USED FOR THIS ANALYSIS 


The program was developed allowing the user to choose from three hull forms to 
run the analysis. The three hull forms are: a basic barge (rectangular hull form), the 
analytical Wigley hull form, and a file containing the data of an existing hull form. 

The barge hull form is used as a demonstration tool for the program. The 
rectangular barge will have the dimensions L (length), B (Beam), T (draft), and H 
(height). The rectangular shape allows the programmer to develop a simple model, 
which can be verified with hand calculations. The user must define all of the necessary 
dimensions at the beginning of the simulation. 

The Wigley hull form used in the development and testing of the program is 
shown in (Figure 2.1). It was chosen due to its ease of analytical representation. The 


offsets of the Wigley hull are described by the following equation. 


ag al a pany 
aE 


Where: 


x = longitudinal distance from midships B = beam (maximum) 
T = draft (maximum) 
L = length between perpendiculars 


y = transverse distance from centerline 


z = height above keel 


This equation for the offsets can be manipulated analytically to derive all of the 
necessary equations to model flooding in this hull. The choice for the rectangular barge 
and Wigley hull forms, besides the obvious advantages of ease of analytical 
representation, offers another feature, which is proven very useful in the parametric 
studies. The barge is a perfectly wall-sided hull form whereas the Wigley hull shown in 
(Figure 2.1) offers a monotonically increasing flare throughout the range of its length. As 
a result it has a significantly higher stiffness in transverse motions. The two hulls can be 
viewed as rather extreme cases, which means that in a qualitative sense, results obtained 


for realistic hull forms are expected to lie between these two cases. 
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Figure 2.1 View of Wigley Hull 


The third hull form will be a more specific user supplied hull. The user must 
supply four files containing the necessary stability parameters to run the simulation. One 
file will contain data for the main hull and the other three files will contain data for the 
three flooded compartments. All four files must be saved in a comma delimited format 
(CSV) and must be ordered such that draft or compartment volume increases until file 
termination. The hull data file must contain six columns of data, where each column 
respectively represents the draft, volume, tons per inch immersion, vertical distance of 
the center of buoyancy above the keel, and the transverse and longitudinal metacentric 
radii. Each of the compartment data files must be formatted with seven columns that 
respectively represent water level, volume, vertical center of gravity, transverse moment 
of inertia, longitudinal moment of inertia, surface area and longitudinal center of gravity. 
This option gives the user the capability to simulate any hull form that has the necessary 
parameters available. For the purpose of program testing and comparison four Matlab M- 
files listed in (Appendix A) were used to generate four files using the analytical Wigley 
equations. A sample of the data generated for the Wigley hull is listed in (Appendix B) 


and a sample of the data for a compartment in the Wigley hull is listed in (Appendix C). 


C. FLOW RATE THROUGH HOLES 


Holes can be modeled as a short-tube orifice with the diameter (hole size) much 
larger than the tube length (hull thickness). This enables the flow through these holes, 
either primary or secondary flooding, to be modeled as turbulent flow through an orifice 


of negligible length. The flow rate then becomes a function of the hole size, the hole 


shape and the pressure difference across the hole. The hole shape is used to derive a 
discharge coefficient Cz. The coefficient, Cz, is supplied by the user at the beginning of 
the simulation. A C, of 0.7 is supplied as the default value based upon predictions for a 
sharp edged hole [Ref. 4]. The equation for flow through the hole can be described by 


the following equation. 


OC, Ale Ne (2) 
Where: 
O = flow rate g = gravitational constant 
A = cross-sectional area of the hole h = the height of fluid on each side of 
Cy = discharge coefficient the hole 


The difference in height of the flooding water on each side of the hole will depend 
on the draft of the ship, trim of the ship (for both exterior water and tank levels) and the 
level of flooding water inside the compartment. In the progressive flooding calculations, 
equation (2) is in reality a set of equations and are implemented as follows: 

e First, the flow rate for each hole in each compartment is calculated using 
equation (2). The square root indicated in the equation, is a signed square root, 
in other words the square root of the absolute value of the head difference Ah 
divided by the sign of Ah. In this way, the value of Q is either positive 
(indicating net inflow into the compartment) or negative (indicating net outflow 


from the compartment) as it should. 


e The head difference Ad depends on the height of the water inside the 
compartment which is calculated by applying the incompressible flow 
continuity equation » Q=A-h,. In this equation, the sum on Q contains all 
inflows and outflows from each compartment including pumping rates, Ac is 
the plan form area of the compartment at the indicated height, and hc is the 
water height inside the compartment. 

e The head difference also depends on the relative position of the damage hole 
with respect to the compartment water level and the seawater level outside the 
hull. The latter depends on the ship’s draft and trim, which are calculated as 
shown in the following section. The ship’s trim also affects the water height 
between adjacent compartments. This trim effect, although relatively small, 


has been incorporated in the simulation program. 


D. NAVAL ARCHITECTURE EQUATIONS 


Regardless of the type of hull-form chosen to simulate there are numerous 
stability parameters that are common to the three hull-forms and these parameters must 
be calculated to simulate the overall effect on stability. Each of the parameters will be 
described below in the order, which they were solved in the model. The computational 


process for each parameter will be explained for each hull-form. 


1. The Distance of the Center of Buoyancy Above the Keel (KB) 


Formally this value can be calculated with the expression 





cian (3) 


where M,, is the total moment of the displaced volume and V is the displaced volume 
[Ref. 5: p.44]. For barge calculations this is always 


— T 
KB =—. (4) 


For the Wigley hull equation (3) must be solved analytically. The value for the total 
moment of the displaced volume must be solved using the expression 


My, = |Zav (5) 


which equates to 





Moc = fof" (7 2) ple rk _ 6 


where x,, Xy are the limits of integration in the longitudinal direction, zy is the limit of 
integration in the vertical direction, Z is height above the keel, L is length between 
perpendiculars and 7 is draft. The displaced volume of the Wigley hull, which will be 


used for numerous calculations equates to the analytical form 


V= [| [dydzde. (7) 


x yg 


Which simplifies to 


ae 4 ees 
vot selon SF-#)) . 


As computed the displaced volume will be in units of ft?. The volume can be converted 


to displaced weight by using the equation 


3 
en = 
35 F 
LT 





(9) 


where A is the displaced weight in long tons. Using the above equations (3), (6) and (8) 


KB simplifies to 











& =H 

2 Sime 

0 = ees (10) 
es 
z i, 


For the user supplied hull form both the value of KB and V must be given as a function 


of draft (7) or simplified as KB(T) and V(T). 


2. Transverse Metacentric Radius (BM) 


Formally this value can be calculated by 


— = 


B Sr (11) 


where /7 is the transverse moment of inertia of the area at the waterline [Ref. 5:p. 84]. 
For the barge the simplifies to 


—— 


& 





. (12) 


where B is the beam and 7 is draft. For the Wigley hull the value /7 is needed. J; can be 


calculated from the relation 


I, = | [x?dA, (13) 


which simplifies to, 


2 


2 3 
I, =B , = 2n) x Ax 42 Ass 4 ey? 
T 5 7 





} (14) 


Where A=4/L’ and xy, x, are the limits of integration in the x (longitudinal) direction. 


With both /7 and V calculated equation (11) can be used. The user supplied hull form 


will contain the relation BM (T). 


3. Tons Per Inch Immersion (7P/) 
This value defines the number of tons of additional weight added that is required to 


submerge the hull one additional inch. This value can be calculated from the expression 


TP] =~, (15) 
420 





where A,, is the area of the waterline [Ref. 5:p. 46]. For the barge the value simplifies to 


* 
TP] =2—.. (16) 


For the Wigley hull A,, must be calculated before 7PJ can be solved. The calculation for 


Ay assumes the analytical form 


B (T-Zy Ae x: 
4,> [fart =a 2) - — [xe-».-4{-4)} (17) 


With A,, calculated equation (15) can be used to solve for 7P/. The user supplied hull 





form will contain the values 7P/(T). 


4. Compartment Volumes (V.) The flooded volumes, of each of the 
compartments, as flooding progresses needs to be computed. For the barge hull 


this computation is simply 


Van), (18) 
where / and b are the length and breath of the compartment and 4 is the height of water in 
the compartment. To compute this for the Wigley Hull equation (7) was used with the 
values of the integration along x equal to the forward and aft bulkheads of the 


compartment. For the user supplied hull type the values of V.(h) must be given. 


5. Parallel Sinkage (PS) 
With the displaced weights of each of the compartments previously calculated the 


parallel sinkage of the hull due to flooding can be found by using the expression 


LA. 
PS =. (19) 
TPI 


This expression can be used for all three hull types using the calculated or supplied 


compartment volumes. 


6. Vertical Height of the Center of Gravity Above the Keel (KG) 
For all hull types an initial KG must be given to start the stability calculations. 


With the initial KG given, the change to KG can be calculated using the following 


formula 


7G = KG" Bola. he.) (20) 


For this calculation the vertical center of gravity for each flooded compartment must be 


calculated. For the barge that value is 


kg. =, (21) 


where / is the height of water in the flooded compartment. For the Wigley hull the center 


of gravity of each compartment can be calculated using the expression 











: [zdv 
Ec a ve > (22) 
which simplifies to 
— 27 ee. 4 1 
kg, = o( 2 - ae \ —X, ~ aa lei eu +}. (23) 


For the user supplied hull form the value kg. (T ) must be given in the data file. Once the 
values of each compartment center of gravity are found equation (20) can be used to find 


the entire hull center of gravity. 


7. Free Surface Correction (FSC) 

When a compartment floods there is a virtual rise in the weight of the liquid in the 
compartment due to the free surface of the liquid. This virtual rise produces a virtual 
moment of free surface equal to the product of the weight of the liquid in the tank and its 
virtual rise. This moment of free surface affects the transverse stability of the ship by 
decreasing the ship's ability to right itself. The free surface correction can be calculated 


as follows, 
GG, = i (24) 


where i; is the transverse moment of inertia of the compartment, 6, is the density of the 


fluid that the ship floats in and 6,1s the density of the fluid in the compartment. The 


16 


transverse moment of inertia must be calculated for each flooded compartment. The 


general form for a rectangular transverse moment of inertia is 


Ib? 


ape (25) 


= 


where / is the length of the compartment and 6 is the beam of the compartment. For the 
barge equation (25) can be used directly to solve for i,. For the Wigley hull equation (14) 
can be used to find the i, of each compartment by changing the values of integration 
along the longitudinal (x) direction to match the values of the forward and aft bulkhead in 
each compartment. The combined free surface correction is the sum of the free surface 


correction terms for each compartment expressed as 


FSC = > gg, , (26) 


where gg, is the free surface correction terms for each compartment. 


8. Metacentric Height (GM ) 
The righting arm of a ship, GZ , is a measure of a ship’s ability to right itself 
when it experiences heel. GZ can be calculated from the formula 
GZ =GM sin d, (27) 
where ¢ is the angle of heel. Equation (27) is not used in the simulation since the model 
does not experience an angle of heel, therefore the metacentric height is used as a direct 


measure of the ship’s ability to right itself. The metacentric height is calculated from the 


following equation for all three hull forms 


GM = KB+ BM — KG-FSC , (28) 


equations (3), (11), (20) and (26) are used for this calculation. The user can predetermine 


a minimum value for GM that will terminate the simulation, signifying a catastrophic 


loss of stability. 


9, Longitudinal Metacentric Radius (BM, ) 

The understanding of the effects of trim and the definitions associated with trim is 
essential, before attempting to determine the changes in draft that results from the 
addition of weight, that occurs during flooding. The longitudinal metacentric radius is 
derived the same as the transverse metacentric radius, by integrating longitudinally to 
obtain the volumes and moments of volumes of the emerged and immersed wedges [Ref. 


5:p. 137]. The formal equation for the longitudinal metacentric radius is 





BM, =—, (29) 


fe 


where /; is the longitudinal moment of inertia of the waterplane. For the Barge hull form 


this calculation takes the form 








if: 
BM, = 30 
LET (30) 
For the Wigley hull form /; must first be found. /; equates to 
1, = [[x?a, (31) 


which simplifies to 


I= fi Copel [ah Ay en } (32) 


Once J; is found equation (29) can be used to solve for BM, . For the user supplied hull 


for the relation BM, (T) must be given. 


10. Longitudinal Free Surface Correction (FSC;) 
The longitudinal free surface correction must be calculated for all flooded 
compartments. Equations (24) and (26) can be used for this calculation if i, for each 


compartment is replaced with the correct longitudinal calculation i. 


11. Longitudinal Metacentric Height (GM, ) 
The longitudinal metacentic height for all three hull forms can be calculated from 


the following expression. 


GM, = KB+ BM, -KG- FSC, (33) 


12. Moment to Change Trim One Inch (M7/) 

The moment to change trim one inch is a convenient quantity to calculate 
because it is independent of trimming moment and it can be precalculated and treated as a 
property of the hull form. The property MT/ can be calculated from the following 


expression for all three hull forms. 


argo 
ye 





(34) 


E. SIMULATION PROGRAM OPERATION 


When first entering the simulation program the main run time screen will appear 


as shown in (Figure 2.2). 
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Figure 2.2 Run Time Environment Screen 


This screen will give the user the option to choose which hull form is desired to run the 
simulation. The status bar provides the user with hints for follow-up actions. From the 
tool bar at the top of the screen one of the three available hull forms must be chosen. 
Some of the menus and buttons on the toolbar are disabled until a hull form is selected. 
When a hull form is chosen the simulation data screen will appear on the as shown in 


(Figure 2.3). 
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Figure 2.3 Simulation Data Screen 


The simulation data screen prompts the user to enter the data necessary to simulate 
the flooding process. The user should select from one of the three tabs across the top of 
the screen. The hull data tab should be used to enter the specific hull and compartment 


dimensions such as length, beam, height, and bulkhead locations. The user must also 


enter the initial KG to start the stability calculations. If a user supplied hull form is 


being simulated the user will be prompted to enter the four files as shown in (Figure2.4). 
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Figure 2.4 Data Files Input Screen 


In the case of the rectangular barge or the Wigley hull, these “open dialog boxes” 
do not appear. When the flooding data tab is selected the user must enter the data 
necessary to simulate the flooding process. The data for all holes is entered on this form 
as shown in (Figure 2.5). All units in this version of the program are in feet. Hole data 
include the height above the keel, hole diameter, its longitudinal position measured 
forward of the aft bulkhead (main compartment only), and discharge coefficient Cg. 
Since roll dynamics are not incorporated in this version of the program, the transverse 
locations of the holes in the secondary compartments are not required. A maximum of 


two main damage holes and three fragmentation holes in each compartment is assumed. 
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Figure 2.5 Flooding Data Input Screen 
When the pumps tab is selected the user can input the data necessary for dewatering as 
shown in (Figure 2.6). 
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Figure 2.6 Pumps Data Input Screen 
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A maximum of two positive displacement pumps can be simulated in each 
compartment. The user must enter the necessary data such as pumping rate in gallons per 
min (GPM) and vertical location. After the simulation data has been entered the user will 
be prompted by the run time environment screen shown in (Figure 2.2) to enter the 
simulation parameters by selecting the simulation parameters tool button at the top of the 
screen. After clicking the simulation parameters tool button the simulation parameters 


screen will appear as shown in (Figure 2.7). 
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Figure 2.7 Simulation Parameters Screen 


The user will have the option of selecting between four tabs across the top of the 
screen. The solver tab will allow the user to enter the parameters used by the Euler fixed 
step differential solver used in the simulation. Euler’s method was chosen for its 
simplicity and because it is accurate for the time scales involved in the problem. The 
duration of the simulation can be entered along with the time step used. A typical time 


step of 0.1 minutes was determined to produce accurate results. The initial conditions 


24 


screen, shown in (Figure 2.8) can be used to regulate an existing water level in the three 
compartments. It should be mentioned that entering an arbitrary water level is not 
allowed in this version of the program and the corresponding text boxes shown in the 
figure do not accept user input. The reasons for this is that a non-zero water level should 
be consistent with the ship’s draft and trim and before the simulation starts, no cross 
check is done by the program. The initial conditions screen can be refreshed if the user 
wishes to continue on a simulation that has ended. In this case the final water levels and 
the ship’s draft and trim are used as initial conditions for the next simulation. The time 
index is also adjusted to account for simulation continuation 1n this case. 
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Figure 2.8 Initial Conditions Screen 


The exit criterion screen, shown in (Figure 2.9), allows the user to enter what 
parameters, if any, may force simulation termination. The user has the option to choose 
between available freeboard and metacentric height. These values are continuously 
monitored during the simulation process. Simulation 1s terminated if any of these two 


minimum values are violated. The pump tab is for future model development 
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incorporating pump automatic control logic and has not yet been incorporated in the 
simulation process. All pumping actions in this version are manually controlled through 


the run time environment screen that is discussed below. 
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Figure 2.9 Exit Criterion Screen 


When the simulation parameters have been the entered the user will be prompted 
to select a file to save the simulation results. This file can be used within the simulation 
program to graph the results or it can be exported in a comma separated values (CSV) 
format, which is compatible with most spreadsheet programs. After the file has been 
selected the main simulation screen will appear as shown in (Figure 2.10). The user must 
click on the run button located at the toolbar to begin the simulation. A progress bar 
depicts schematically simulation time progress. Also shown are the water levels in the 
tanks, and the ship waterline. Although these are depicted as square, the correct 
geometry as entered is used in the calculations. Damage and fragmentation hole 


locations and heights are also shown. The user may conduct real time damage control 
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operations using the hole clogging and pumps tabs (Figure 2.11) at the bottom of the 
screen. Hole diameters as well as pump heights and settings can be continuously 


adjusted and observed during simulation. 
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Figure 2.10 Main Simulation Screen 
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Figure 2.11 Pumps Manipulation Screen 


The user can monitor progress of the simulation by watching various continuous 


updates on the screen including the water level in each compartment, the simulation time, 


draft forward and aft, available freeboard and available GM . After successful program 
termination the user will have the options to continue the simulation, export the data in 
comma separated format, or graph the results. All main results can be also presented in 
graphical forms using the program’s built-in graphical capabilities. A sample graph of 
the water depth of all compartments at simulation termination is shown in (Figure 2.11). 
These results (both in raw data or graph) can be copied into the Windows clipboard and 


pasted into any other application. 
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Ul MODEL VALIDATION 


When conducting model validation various calculations of the model were tested 
separately. First the validation of the flow rates entering the main compartment were 
tested and compared to U. S. Navy standards. Next the effects of flooding on each of the 
three hull forms were observed by using various scenarios to demonstrate the progressive 
flooding and damage control process. The scenarios were used to both validate the 
model and analyze the hull forms’ transient and steady state characteristics of progressive 


flooding. 
A. FLOW RATES VALIDATION 


For validation of the flow rates of the flooding in the main compartment the U.S. 
Naval Ship’s Technical Manual (NSTM) [Ref. 6] was used. Table 1 below shows both 
the NSTM predictions for flow rate and the model observed values with a Cz of 1.0 used 


as an input parameter and a hole diameter of six inches. 


— eee | 
NSTM | Q(gal/Min) | 1000 | 1414 [1732 | 2000 | 2236 | 2449 2828 | 3000 
ee | | 


Table 1. Comparison of Model and NSTM flow rates 









3) 


As can be seen from Table 1 the model accurately calculates the flow rate through hole 


with a percent difference from the NSTM values of less than 0.25%. 


B. FLOODING ANALYSIS 


To demonstrate the ability of the model to accurately model progressive flooding 
and to analyze some of the transient and steady state characteristics of progressive 
flooding various scenarios were developed and ran. The scenarios will be run on the 
three different hull forms to compare the similarities and differences. 

The barge hull form was first tested as scenario 1. A 120-minute scenario was 
run using all default parameters. No damage control or dewatering equipment was used 
because only correct model operation was being evaluated. The main simulation screen 
at scenario termination is shown in (Figure 3.1). As can be seen from the figure the main 
compartment had a high flooding level and the secondary compartments also experienced 


progressive flooding as expected. 
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Figure 3.1 Scenario 1 Barge Hull Simulation Screen 


DZ 


Other results must be graphed to further evaluate the results. The flow rates for 
the main compartment are shown in (Figure 3.2). The graph demonstrates the expected 
results. Both holes have the correct initial flow rate, which was previously proven, and 


the flow rate continues to increase due to hull sinkage until the hole becomes submerged. 
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Figure 3.2 Scenario 1 Main Flow Rates vs. Time 


As can be seen from (Figure 3.3) the water level in the main compartment 
increased sharply at the beginning of the simulation. The rate of increase is nearly linear, 
as expected, because of the nearly constant flow rates until the hole becomes submerged 


shown in (Figure 3.2). The nearly linear region of the graph ends and the slope begins to 
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decrease when the main holes become submerged and progressive flooding takes effect. 
The water levels in the fore and aft compartments indicate the occurrence of progressive 
flooding. Progressive flooding can not begin until the level of water in the main 
compartment reaches the height of the secondary fragmentation holes. This occurs at 
approximately the 20-30 minute range, when the water level in the main compartment 
exceeds five feet, which is the height of the lowest fragmentation hole. The water level 
in the forward compartment increases at a faster rate than the aft compartment. This is 
because the forward compartment contains larger fragmentation holes than the aft 
compartment. 

Shown in (Figure 3.4) are the forward and aft drafts of the barge hull as the 
simulation progressed. The forward draft increased at a faster rate than the aft indicating 
a trim by the bow condition. This was expected because the compartments are located 
forward of midships. 


The transverse metacentric height was evaluated in (Figure3.5). The effects of the 


free surface correction to GM can be seen by the sharp decrease. Free surface correction 
has a drastic correction on the barge hull form because the rectangular surface area of the 
flooded compartment starts at a large value and remains constant. The Wigley hull form, 
which is more representative of actual hulls, will not show this trend because the surface 
area of the flooded compartment starts at a low value and increases as the level of 
flooding increases. 

Effects of flooding that are common to all hull forms, such as: initial increase in 


hull flow rate due to increasing depth of the hole; decrease in flow rate due to the 
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submergence of the hole; etc., will not be readdressed in following scenarios unless the 


effect is specific to the type of hull being addressed. 
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Figure 3.3 Scenario 1 Water Levels vs. Time 
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Figure 3.4 Scenario | Hull Draft 
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Figure 3.5 Scenario 1 GM vs. Time 


Scenario 2 was ran to demonstrate flooding effects on the Wigley hull form. The 
scenario ran with all default parameters for 120 minutes. No dewatering equipment was 
used. All useful parameters, such as, flooding rates, final draft, final flooding levels, and 
metacentric height will be shown so that a comparison with the user supplied data hull 


form can be made in scenario 3. The main simulation screen at scenario termination is 


oy 


shown in (Figure 3.6). As can be seen from the figure the main compartment had a high 
flooding level and the secondary compartments experienced progressive flooding as 


previously demonstrated by the barge hull form. 
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Figure 3.6 Scenario 2 Wigley Hull Simulation Screen 


As shown in (Figure 3.7) the water levei in the main compartment quickly reach 
an equilibrium level similar to the barge followed by progressive flooding to the fore and 
aft compartments. There are two notable differences between the flooding levels in the 
barge and the Wigley hull form. First, there is a high initial flooding level in the main 
compartment. This occurs because the Euler ODE solver calculates the level of water in 


each compartment from the following equation 
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Level (ft) 


= Q 
=n, +a 2), (35) 


where / is water height, df is the time step, Q 1s the flow rate, and A is the compartment 
surface area. The surface area of the Wigley hull is zero when the water level is zero and 
increases as water level increases. This makes the initial water level high but has no 
effect on the proper program operation. This occurrence can be corrected by using a 
smaller time step. Second, the water levels in the main and fore compartments reach 
equilibrium unlike the barge hull form. This can also be attributed to the small initial 


compartment surface areas. 


Compartment Water Level vs. Time 
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Figure 3.7 Scenario 2 Water Levels vs. Time 
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GM (ft) 


The GM is shown in (Figure 3.8). The plot shows a steady decrease unlike the 
sharp decrease experienced in scenario 1. The gradual decrease is due to the free surface 
correction. The free surface correction did not have as large an impact as scenario 1 
because, unlike the barge, the surface area in the Wigley hull was not large at low water 


levels. 


GM vs. Time 
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Figure 3.8 Scenario 2 GM vs. Time 


Scenario 3 was run on the user supplied hull form using the data generated by the 
Matlab files in (Appendix A). The purpose of this scenario was to test the accuracy of 


the calculations using a user supplied hull form and compare the results to scenario 2. As 
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expected scenario 3 showed the same trends as previously discussed in scenario 2. Table 


2 lists the tabulated results of scenarios 2 and 3, along with the percent difference. 


FWD Draft | AFT Draft Main Fwd Aft 
Comp Comp Comp 
Level Level Level 
23.24 21.94 DS] 22.65 
0.43% 0.64% 0.09% 0.00% 3.26% 1.54% 


Table 2. Comparison of Scenario | and Scenario 2 Results 












Relatively low percent differences were encountered in all the tabulated results with the 
highest error being the water level in the aft compartment. 

Scenario 4 was ran to demonstrate the effects of hole plugging and dewatering 
while running a simulation on the Wigley hull form. A Wigley hull with default 
dimensions was given standard battle damage as experienced in scenario 2. All pumps 
were activated, but were not started until flooding progressed. Two 3,000 GPM pumps 
were in the main compartment and two 2,000 GPM pumps were placed in each the 
forward and aft compartment. The water heights in each compartment, along with pump 


activation and hole clogging times, are shown in (Figure 3.9). 


4] 


Level (ft) 


Compartment Water Levels vs. Time 


Hole 1 repaired 
at 41 min 


Both pumps 
started 3000 GPM 
each at 10 min 


Fore water Both pumps started 2000 
level GPM each at 45 min 





Time (min) 


Figure 3.9 Scenario 4 Water Height vs. Time 


The slopes of the water level lines in each compartment decreases when each pump is 
activated as expected. The times of activation are indicated on the graph. The water 
level in the main compartment begins to decrease when the hole is plugged. It should be 
noted that the initial water level in the main compartment was lower than scenario 2 
because a smaller time step was used. The flow rates for holes in the main compartment 


are shown in (Figure3.10). The flow rates are shown to demonstrate the decrease when 
hole 1 is plugged at 41 min. The effects on GM are shown in (Figure 3.11). 


GM decreases as flooding progresses then begins to increase while dewatering is being 


conducted as expected. 
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Figure 3.10 Scenario 4 Flow Rates vs. Time 
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Figure 3.11 Scenario 4 GM vs. Time 
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IV. CONCLUCIONS 


This thesis has successfully developed a progressive flooding model using the 
Visual Basic programming language. This model can be used as a design tool to aid the 
engineer in stability verifications following battle damage with progressive flooding. The 
engineer can conduct real time scenarios using the installed dewatering systems or 
damage control hole plugging techniques to better evaluate any hull form. 

Based upon the results discussed in the preceding scenarios the following 

conclusions are drawn: 

1. This model can be used as a design tool to evaluate progressive flooding. The 
engineer can conduct real time scenarios to verify ship stability and 
survivability. 

2. The user can adapt the model to any hull form and conduct scenarios on any 
compartment in the hull. 

3. Vital Naval architectural parameters, such as free surface effects and internal 
tank trim are included in the calculations to give more realistic results. 

4. The user can use multiple pumping combinations and repair or plug holes to 


aid in the damage control process. 
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V. RECOMMENDATIONS 


Because this model was developed from ground zero and due to time scope 
limitations some features of this design tool are still undeveloped. Some of the areas for 
future development are: 
1. Conduct further testing with an actual data supplied hull form with known 
stability parameters and compare the results. 

2. Expand the dewatering feature of the model to include logic controlled 
pumping. 

2. Expand the model to include longitudinal bulkheads so that transverse 


stability criterion, such as list can be included in the simulation. 
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APPENDIX A. 


HULL AND COMPARTMENT DATA GENERATION PROGRAMS 


ia SpLcer 


aoc oO 


This Matlab program computes various hull characteristics for the 
Wigley hull 
% form and saves them to a file named 'hulldata’. 


ao 


clear all 


Z1=0.1; $ Temp initial draft 

Zh=40; $ Temp final Draft 

B=80; $ Breadth 

T=40; % Draft 

X1=-200; $ End of ship (length) 
Xh=200; $ begging of ship (length) 
L=400; $ Length of ship 

KG=25; 

delta=.1; 


hulldata = zeros(200, 6); 


S$Hulldata 1 = draft 
$Hulldata 2 = volume 
S$Hulldata 3 = TPI 
S$Hulldata 4 = KB 
S$Hulidata 5 = BM 
$Hulldata 6 = BML 

A = 4/L%2; 


MEMS siemens 3) toe Ae (xh) 5) fo) (A 3* (xt) oid aS 
GA eal aor reo a eet IES) 0/7 a) (ASS * x 7 7d): 


i=1; 
for 2-21: deltacan 


hudidata (1,1) >= 22. 
hudi data (i) 7) ee— eee ei eee aes ae) (XX I) = 


G47 3) (kno 3) 7 bee) =e sy 2 ey 
Rieliieiat a(S eee 2 ee nee) 2 a i ay) X37 2 
(le 3/712) )))/420; 
no lidatats , 4) = (e427 1) a 2a47 (49992) (2°27 H=2°3/ (3*T*2) ); 
neielieata (io) = (1/3) * Bos a(t 2372/7 2)°3~ (Temp) )/hulidata(i; 2) ; 
huiidatati, 6) = te Gl E-z) 427 172) pet 2) *373—(-L/2)*3/3- 
(Ay (52) ) * ( (72) 2 5—(—-o, 2) 5) ) ol data (i; Z) ; 
ee 
end 


save hulldata hulldata -ASCII; 
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oe 


Tim Spicer 


This program computes various characteristics for the fwd compartment 
i 
the Wigley hull form. 


ce O of oP 


clear all 


Zl Olle: 6 Temp initial draft 

Zh=40; @ Temp final Draft 

B=80; 6 Breadth 

T=40; 6 Draft 

X1=25; % End of compartment (length) 
Xh=40; ® begging of compartment (length) 
L=400; @ Length of ship 

KG=25; 

delta=.1; 


fwddata = zeros(200,7); 


$fwddata 1 = draft 
$fwddata 2 = volume 
$fwddata 3 = zcoord 
$fwddata 4 = It 
$fwddata 5 = Il 
$fwddata 6 = Area 
$fwddata 7 = xcoord 
A = 4/L%2; 


Tempe exe (iene Sar (3 Ae xh) 5) or ( (AOS (xh) OF) / 7) x2 j= 
CPS Kae et es rk oe el) 7/4 (AS (IO ef cs 


iLf=. i 


for z=Zl:delta:Zh 

fwddatalien) =| z- 

fwddatat,2)e— 8 G2 2e tt) (zee rd 2) x) — 
(047 3) en 307 Bo 2) ea 2) 

fuddeta (iss \e—9B* (2. zs (oe peed 7 2) <i h—— 
GAV(37E C2) ) = (xno 3-xd 3) (fwoeatait, 2) 

fwddata(i,4) = ((1/3)7*B*3* (1-(1-z) *2/1-2) *3*% {xh 
(R= 3) +( (S*AC2* in) 75) Jo) Ass ti ee) — 
(Ae x13) + (1(3* ACoA) 5) 7 5) = (AOS Oe ae ee 

fwcdatai1, 5) = (Bs Wai (E-2)e2/ 2°72) 1 (xh 3-3/3 
(AV 5-2) (xXho>5=x1°5)) )e 

Ewadaea(l, 6) = 224872)" (tT TP 2 eat xX l= (4/3) * ns Le 
(cls een): 

fwodataviy)) = (U(x 2—iee 2 2 rn 4 x (2 ee xh 3 = 
Kies) 7 oes 2.) i; 

L=it+l1; 


et le 
end 


save fwddata fwddata -ASCII; 
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ST imsopLcer 

$ 

% This program computes various characteristics for the main 
compartment of 

$ the Wigley hull form. 


clear all 


Z21=0.1; $ Temp initial draft 

Zh=40; $ Temp final Draft 

B=80; $ Breadth 

T=40; %$ Draft 

X1=-10; 6 End of compartment (length) 
MR=Z25% $ begging of compartment (length) 
L=400; $ Length of ship 

KG=25; 

delta=.1; 


mMaindata = zeros(200,7); 


$maindata 1 = draft 
$maindata 2 = volume 
$maindata 3 = zcoord 
$maindata 4 = It 
$maindata 5 = Il 
$maindata 6 = Area 
$maindata 7 = xcoord 
A = 4/L%2; 


Temp = Xh-(A*Xh*3) + ( (3*A*2* (Xh) *5) /5) - ( (A*3* (Xh) *7) /7) - ((X1) - 
(dhe Shc (3 eh 2 sl a) />) = (aes (x) Ae); 


i = 1; 


for z=Zl:delta:Zh 

maindata(i,l) = 2; 

Mamncdata (a2) — Bat ( (2-2/1) = (2-3) 7 STZ) ) ) * Ci xehex — 
i dmaneeny (Miro 87 2 or XS) fe) i) 

Maindata(i,3) =] B* (( ((2*273) 7 (3°77) )=(2%4/7 (47172) ) )* (xh-X1- 
(27 }S41.°2)) * (ene s—ml 3) )/ (maandaca(1, 2) 

maindata(i,4) = ((1/3) *B*3* (1-(T-z) *2/T*2) *3* (Xh- 
(Anco ( (3 *AC2 exh) 75) 75) = (AOS (XH) COM 7) = (XL) 
(Ax) ~Coet ( (S*AC2* (XL) 5) (So (As ial) A777) 

Maindata(i,5) = Bet ( (t=2) 427772) ) kn 3/ 3-—xX1 37 3- 
(7 does 2) )* (Xe oaee oan. 

meneatatl 6) = mwaaiB/ 2) “a Uh 2) 7 2 ( Xe X42 78 x3 / Lo 2 
elas eo} ).) ; 

Maindata(i,7) = (((Xh*2-X1%2) /2)- (Xh*4-X1%4) / (L*%2) ) /( (Xh-X1)- 
mines hes) (3. (ee) ) 

i=it+l; 


i=itl; 
end 


Save maindata maindata -ASCII; 
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Pim “Sp leer 


This program computes various characteristics for the aft compartment 


Fh 


0° O of oP? of 


the Wigley hull form. 


clear all 


Z1=0.1; %¢ Temp initial draft 

Zh=40; $ Temp final Draft 

B=80; 6 Breadth 

T=40; $ Draft 

X1l=-35; * End of compartment (length) 
Xh=-10; % begging of compartment (length) 
L=400; $ Length of ship 

KG=25; 

delta=.1; 


aftdata = zeros(200,/7); 


saftdata 1 = draft 
$aftdata 2 = volume 
$aftdata 3 = zcoord 
S$aftdata 4 = It 

 $aftdata 5 = I1 
S$aftdata 6 = Area 
Zaftdata 7 = xcoord 
A = 4/L%2; 


Temp = xh-(A* Xho sie tones (XN) oD) 75) 1 (A 3* (Xn) 77) 7/7) —( (XT) — 
el) (2 Re a as) 7) a (AOS (XE) a7 77); 


1 = i; 


for z=Zl:delta:Zh 
aftdata(i,1) = 2; 
aAttdataliec) = BF (((2°2)7 (lz 3 ede x) = 
(4/3)? CC (4he8) /L° 2) Xt s hy 2 ea 
aftdatali;3) = Bat(24 27 2)0 sr 2047 (a2 )).) © (xn 
(403 * G42) ) xno 3=xX1 3) (attcatace ae 
attdata(i,4) = ((1/3) "83s (l—(T-2) 2 xh 
(A*Xh*3) + ((3*A%2* (Xh) *5) /5) -( (A*®3* (Xh) *7) /7) -((X1)- 
(Ae Mie) (3 * A 2* (XE) = WA 7) 
atedatatip5s) = (B*(i—((T-2) 7270-2) )* (4he37s_-x1°3/3- 
(AGS ie 2 ye (Xho 5 =X2 5) De 
attdaga(i,6) = 2*(B72)* (lor) T) 720-1 (47 3) (x8 
(es Eazy 
aftdata(i,7) 
RAS tse (eZ), 
1=itl; 
end 


(( (Xno2-K192)7 2) (a4 — x1 a 2 Gn 


save aftdata aftdata -ASCII; 
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APPENDIX B. 


SAMPLE HULL DATA 


1.0000000e-001 5.3288889e+000 2.5365079e-001 6.6659716e-002 7.2929689e-001 1.5993328e+005 

5.0666667e-001 1.3330551e-001 1.4543370e+000 7.9933222e+004 
7.5904762e-001 1.9993734e-001 2.1751303e+000 5.3266500e+004 
2.8916866e+000 3.9933110e+004 
3.6040158e+000 3.1933054e+004 
4.3121279e+000 2.6599665e+004 
5.0160328e+000 2.2790085e+004 
5.7157402e+000 1.9932886e+004 
6.4112602e+000 1.7710607e+004 


2.0000000e-001 
3.0000000e-001 
4.0000000e-001 
5.0000000e-001 
6.0000000e-001 
7.0000000e-001 
8.0000000e-001 
9.0000000e-001 
1.0000000e+000 
1.1000000e+000 
1.2000000e+000 
1.3000000e+000 
1.4000000e+000 
1.5000000e+000 
1.6000000e+000 
1.7000000e+000 
1.8000000e+000 
1.9000000e+000 
2.0000000e+000 
2.1000000e+000 
2.2000000e+000 
2.3000000e+000 
2.4000000e+000 
2.5000000e+000 
2.6000000e+000 
2.7000000e+000 
2.8000000e+000 
2.9000000e+000 
3.0000000e+000 
3.1000000e+000 
3.2000000e+000 
3.3000000e+000 
3.4000000e+000 
3.5000000e+000 
3.6000000e+000 
3.7000000e+000 
3.8000000e+000 
3.9000000e+000 
4.0000000e+000 
4.1000000e+000 
4.2000000e+000 
4.3000000e+000 
4.4000000e+000 
4.5000000e+000 


2.1297778e+001 
4.7880000e+001 
8.5048889e+001 
1.3277778e+002 
1.9104000e+002 
2.5980889e+002 
3.3905778e+002 
4.2876000e+002 
5.2888889e+002 
6.3941778e+002 
7.6032000e+002 
8.9156889e+002 
1.0331378e+003 
1.1850000e+003 
1.3471289e+003 
1.5194978e+003 
1.7020800e+003 
1.8948489e+003 
2.0977778e+003 
2.3108400e+003 
2.5340089e+003 
2./672578e+003 
3.0105600e+003 
3.2638889e+003 
3.5272178e+003 
3.8005200e+003 
4.0837689e+003 
4.3769378e+003 
4.6800000e+003 
4.9929289e+003 
5.3156978e+003 
5.6482800e+003 
5.9906489e+003 
6.3427778e+003 
6.7046400e+003 
7.0762089e+003 
7.4574578e+003 
7.8483600e+003 
8.2488889e+003 
8.6590178e+003 
9.0787200e+003 
9.5079689e+003 
9.9467378e+003 
1.0395000e+004 


1.0107937e+000 
1.2619048e+000 
1.5123810e+000 


2.665551 8e-001 
3.3315900e-001 
3.9974874e-001 


1.7622222e+000 4.6632439e-001 


2.0114286e+000 


5.3288591e-001 


2.2600000e+000 5.9943325e-001 


2.5079365e+000 
2.7552381e+000 
3.0019048e+000 
3.2479365e+000 
3.4933333e+000 
3.7380952e+000 
3.9822222e+000 
4.2257143e+000 
4.4685714e+000 
4.7107937e+000 
4.9523810e+000 
5.1933333e+000 
5.4336508e+000 
5.6733333e+000 
5.9123810e+000 
6.1507937e+000 
6.3885714e+000 
6.6257143e+000 
6.8622222e+000 
7.0980952e+000 
7.3333333e+000 
7.5679365e+000 
7.8019048e+000 
8.035238 1e+000 
8.2679365e+000 
8.5000000e+000 
8.7314286e+000 
8.9622222e+000 
9.1923810e+000 
9.4219048e+000 
9.6507937e+000 
9.8790476e+000 
1.0106667e+001 

1.0333651e+001 

1.0560000e+001 

1.0785714e+001 
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6.6596639e-001 

7.3248528e-001 

7.9898990e-001 

8.6548020e-001 

9.3195616e-001 

9.9841772e-001 

1.0648649e+000 
1.1312975e+000 
1.1977157e+000 
1.2641194e+000 
1.3305085e+000 
1.3968830e+000 
1.4632428e+000 
1.5295879e+000 
1.5959184e+000 
1.6622340e+000 
1.7285349e+000 
1.7948210e+000 
1.8610922e+000 
1.9273484e+000 
1.9935897e+000 
2.0598161e+000 
2.1260274e+000 
2.1922237e+000 
2.2584048e+000 
2.3245708e+000 
2.3907216e+000 
2.4568573e+000 
2.5229776e+000 
2.5890827e+000 
2.6551724e+000 
2.7212468e+000 
2.7873057e+000 
2.8533492e+000 
2.9193772e+000 


7.1026026e+000 
7.7897774e+000 
8.4727943e+000 
9.1516634e+000 
9.8263944e+000 
1.0496997e+001 
1.1163482e+001 
1.1825858e+001 
1.2484136e+001 
1.3138325e+001 
1.3788436e+001 
1.4434478e+001 
1.5076460e+001 
1.5714394e+001 
1.6348288e+001 
1.6978153e+001 
1.7603999e+001 
1.8225835e+001 
1.8843672e+001 
1.9457518e+001 
2.0067385e+001 
2.0673281e+001 
2.127521 7e+001 
2.1873203e+001 
2.2467248e+001 
2.3057363e+001 
2.3643556e+001 
2.4225839e+001 
2.4804221e+001 
2.5378711e+001 
2.5949320e+001 
2.6516058e+001 
2.7078933e+001 
2.7637957e+001 
2.8193139e+001 


1.5932773e+004 
1.4478171e+004 
1.3265993e+004 
1.2240296e+004 
1.1361118e+004 
1.0599156e+004 
9.9324324e+003 
9.3441400e+003 
8.8212070e+003 
8.3533134e+003 
7.9322034e+003 
7.5511935e+003 
7.2048156e+003 
6.8885523e+003 
6.5986395e+003 
6.3319149e+003 
6.0857031e+003 
5.8577247e+003 
5.6460263e+003 
5.4489237e+003 
5.2649573e+003 
5.0928558e+003 
4.9315068e+003 
4.7799330e+003 
4.6372717e+003 
4.5027590e+003 
4.3757159e+003 
4.2555367e+003 
4.1416795e+003 
4.0336580e+003 
3.9310345e+003 
3.8334140e+003 
3.7404392e+003 
3.6517859e+003 
3.5671595e+003 


2.9853896e+000 2.8744489e+001 3.4862915e+003 
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APPENDIX C. 


SAMPLE COMPARTMENT DATA 


1.0000000e-001 4.9260998e-001 3.3329858e-002 5.0950090e-001 5.4716415e+003 9.8480911e+000 -2.2440577e+001 


2.0000000e-001 
3.0000000e-001 
4.0000000e-001 
5.0000000e-001 
6.0000000e-001 
7.0000000e-001 
8.0000000e-001 
9.0000000e-001 
1.0000000e+000 
1.1000000e+000 
1.2000000e+000 
1.3000000e+000 
1.4000000e+000 
1.5000000e+000 
1.6000000e+000 
1.7000000e+000 
1.8000000e+000 
1.9000000e+000 
2.0000000e+000 
2.1000000e+000 
2.2000000e+000 
2.3000000e+000 
2.4000000e+000 
2.5000000e+000 
2.6000000e+000 
2.7000000e+000 
2.8000000e+000 
2.9000000e+000 
3.0000000e+000 
3.1000000e+000 
3.2000000e+000 
3.3000000e+000 
3.4000000e+000 
3.5000000e+000 
3.6000000e+000 
3.7000000e+000 
3.8000000e+000 
3.9000000e+000 
4.0000000e+000 
4.1000000e+000 
4.2000000e+000 
4.3000000e+000 
4.4000000e+000 
4.5000000e+000 


1.9687965e+000 
4.4260945e+000 
7.8620389e+000 
1.2274164e+001 

1.7660006e+001 

2.4017099e+001 

3.1342978e+001 

3.9635177e+001 

4.8891233e+001 
5.9108679e+001 
7.0285050e+001 
8.2417882e+001 
9.5504709e+001 
1.0954307e+002 
1.2453049e+002 
1.4046451e+002 
1.5734267e+002 
1.7516250e+002 
1.9392153e+002 
2.1361730e+002 
2.3424735e+002 
2.5580920e+002 
2.7830040e+002 
3.0171848e+002 
3.2606097e+002 
3.5132541e+002 
3.7750934e+002 
4.0461028e+002 
4.3262578e+002 
4.6155337e+002 
4.9139058e+002 
5.2213495e+002 
5.5378401e+002 
5.8633530e+002 
6.1978635e+002 
6.5413470e+002 
6.8937788e+002 
7.2551344e+002 
7.6253889e+002 
8.0045178e+002 
8.3924964e+002 
8.7893001e+002 
9.1949043e+002 
9.6092842e+002 


6.6652755e-002 
9.9968672e-002 
1.3327759e-001 
1.6657950e-001 
1.9987437e-001 
2.3316220e-001 
2.6644295e-001 
2.9971662e-001 
3.329831 9e-001 
3.6624264e-001 
3.9949495e-001 
4.3274010e-001 
4.6597808e-001 
4.9920886e-001 
5.3243243e-001 
5.6564877e-001 
5.9885787e-001 
6.3205970e-001 
6.6525424e-001 
6.9844148e-001 
7.3162139e-001 
7.6479397e-001 
7.979591 8e-001 
8.3111702e-001 
8.6426746e-001 
8.9741049e-001 
9.3054608e-001 
9.6367421e-001 
9.9679487e-001 
1.0299080e+000 
1.0630137e+000 
1.0961118e+000 
1.1292024e+000 
1.1622854e+000 
1.1953608e+000 
1.2284286e+000 
1.2614888e+000 
1.2945413e+000 
1.3275862e+000 
1.3606234e+000 
1.3936528e+000 
1.4266746e+000 
1.4596886e+000 
1.4926948e+000 


4.0607222e+000 
1.3653480e+001 
3.2242135e+001 
6.2735883e+001 
1.0799903e+002 
1.7085130e+002 
2.5406809e+002 
3.6038076e+002 
4.9247693e+002 
6.5300075e+002 
8.4455319e+002 
1.0696923e+003 
1.3309335e+003 
1.6307499e+003 
1.9715725e+003 
2.3557905e+003 
2.7857516e+003 
3.2637623e+003 
3.7920880e+003 
4.3729536e+003 
5.0085435e+003 
5.7010018e+003 
6.4524330e+003 
7.2649018e+003 
8.1404336e+003 
9.0810149e+003 
1.0088593e+004 
1.1165078e+004 
1.2312339e+004 
1.3532210e+004 
1.4826485e+004 
1.6196923e+004 
1.7645244e+004 
1.9173132e+004 
2.0782232e+004 
2.2474157e+004 
2.4250479e+004 
2.6112738e+004 
2.8062436e+004 
3.0101041e+004 
3.2229985e+004 
3.4450666e+004 
3.6764446e+004 
3.9172655e+004 


oe 


1.0929587e+004 
1.6373836e+004 
2.1804389e+004 
2.7221245e+004 
3.2624405e+004 
3.8013870e+004 
4.3389637e+004 
4.8751709e+004 
5.4100085e+004 
5.9434764e+004 
6.4755747e+004 
7.0063034e+004 
7.5356624e+004 
8.0636519e+004 
8.5902717e+004 
9.1155219e+004 
9.6394024e+004 
1.0161913e+005 
1.0683055e+005 
1.1202826e+005 
1.1721228e+005 
1.2238261e+005 
1.2753924e+005 
1.3268217e+005 
1.3781141e+005 
1.4292695e+005 
1.4802879e+005 
1.5311694e+005 
1.5819139e+005 


1.9671531e+001 -2.2440577e+001 
2.9470320e+001 -2.2440577e+001 
3.9244458e+001 -2.2440577e+001 
4.8993945e+001 -2.2440577e+001 
5.8718781e+001 -2.2440577e+001 
6.8418966e+001 -2.2440577e+001 
7.8094500e+001 -2.2440577e+001 
8.7745383e+001 -2.2440577e+001 
9.7371615e+001 -2.2440577e+001 
1.0697320e+002 -2.2440577e+001 
1.1655013e+002 -2.2440577e+001 
1.2610240e+002 -2.2440577e+001 
1.3563003e+002 -2.2440577e+001 
1.4513301e+002 -2.2440577e+001 
1.5461133e+002 -2.2440577e+001 
1.6406501e+002 -2.2440577e+001 
1.7349403e+002 -2.2440577e+001 
1.8289840e+002 -2.2440577e+001 
1.9227813e+002 -2.2440577e+001 
2.0163320e+002 -2.2440577e+001 
2.1096361e+002 -2.2440577e+001 
2.2026938e+002 -2.2440577e+001 
2.2955050e+002 -2.2440577e+001 
2.3880697e+002 -2.2440577e+001 
2.4803878e+002 -2.2440577e+001 
2.5724595e+002 -2.2440577e+001 
2.6642846e+002 -2.2440577e+001 
2.7558632e+002 -2.2440577e+001 
2.8471953e+002 -2.2440577e+001 


1.6325214e+005 2.9382809e+002 -2.2440577e+001 
1.6829920e+005 3.0291200e+002 -2.2440577e+001 
1.7333256e+005 3.1197126e+002 -2.2440577e+001 
1.7835223e+005 3.2100586e+002 -2.2440577e+001 
1.8335820e+005 3.3001582e+002 -2.2440577e+001 
1.8835047e+005 3.3900113e+002 -2.2440577e+001 
1.9332905e+005 3.4796178e+002 -2.2440577e+001 
1.9829393e+005 3.5689778e+002 -2.2440577e+001 
2.0324512e+005 3.6580913e+002 -2.2440577e+001 
2.0818260e+005 3.7469583e+002 -2.2440577e+001 
2.1310640e+005 3.8355788e+002 -2.2440577e+001 
2.1801649e+005 3.9239528e+002 -2.2440577e+001 
2.2291289e+005 4.0120803e+002 -2.2440577e+001 
2.2779560e+005 4.0999612e+002 -2.2440577e+001 
2.3266460e+005 4.1875957e+002 -2.2440577e+001 
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